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ABSTRACT 



A technique of statistically adjusting dynamical fore- 
casts of tropical cyclone motion was tested. All tests were 
performed with operationally-analyzed data from the U.S. Navy 
Fleet Numerical Weather Central (FNWC). Three sets of regres- 
sion equations were developed to modify forecasts of typhoon 
tracks. The first set of equations was based only on forward 
integration of the FNWC Tropical Cyclone Model (TCM) for 28 
cases in 1975-76. An independent sample of cases from 1977- 
78 indicated that the first equation set was based on too 
small a sample size, especially considering the anomalous 
nature of the 1975-75 storm tracks. A second equation set 
based only on forward integration of the TCM was derived from 
61 storm track forecasts from 1975-78. Results from the ex- 
periments with these equations indicate that systematic data 
and model errors can be used to statistically adjust forecast 
storm tracks. The second equation set based on forward inte- 
gration showed improvement over the unmodified model predic- 
tions at all forecast times. A third equation set based on 
forward and backward integration of the TCM explained the 
greatest amount of variance of all the equation sets. In a 
dependent test of these equations using 31 of the 1977-78 
cases, the U.S. Navy 7th Fleet error goal of 100 and 150 
nautical miles at 43 and 72 hours was nearly met. 
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I. INTRODUCTION 



The tropical storm is one of nature's most destructive 
phenomena, making typhoon p re diction a major concern world- 
wide. The need to forecast storm conditions has produced a 
variety of prediction methods. Subjective methods, includ- 
ing persistence, have generally produced fairly good short- 
range forecasts. An especially difficult problem for any 
subjective or objective prognostic scheme is the occurrence 
of storm recurvature. This problem was significant in the 
western North Pacific Ocean during the 1975 typhoon season, 
and to a lesser extent during the 1976 typhoon season. An 
unusual number of storms tracked northward, while others re- 
curved. As a result, the forecast error statistics were 
higher than normal (Annual Typhoon Report, 1975). The 
United States Seventh Fleet Commander, noting the necessity 
to improve the forecast errors, has levied a requirement for 
the Joint Typhoon Warning Center (JTWC) to achieve maximum 
forecast errors of 50, 100 and 150 nautical miles for 24, 

48 and 72 hours, respectively. In comparison, the 1977 JTWC 
average position forecast errors for tropical cyclones at 24- 
48-, and 72-hours were 140, 266 and 390 nautical miles, re- 
spectively (Annual Typhoon Report, 1977). The long-term 
error trend, as indicated by the five-year mean error, has 
been increasing since 1972. In 1977 the operational TCM 
produced mean vector errors of 138, 262 and 450 nautical 
miles for all tropical cyclones (Annual Typhoon Report, 1977) 
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Although progress has been made in the numerical simula- 
tion of tropical storm characteristics such as intensity and 
spiral bands (Anthes et al, 1974), the major emphasis has 
focused on storm motion (Hovermale et al, 1976; Ley, 1975). 
Multi-level, nested-grid models are being developed in an 
attempt to improve forecast positions over those of analog, 
persistence and statistical techniques. In spite of the 
sophistication of some numerical models, undesirable conse- 
quences may result from initialization with poor or limited 
data. One test of the dynamical models is the accuracy of 
the initial storm track. If the short-term forecasts are 
inaccurate, a dynamic model cannot be expected to produce 
accurate, extended forecasts. Some dynamic models (Hovermale 
et al, 1976; Ley, 1975; Hodur and Burk, 1978) predict tracks 
with systematic bias relative to the actual track and com- 
monly predict motion which is too slow. A likely source 
of error seems to be in specification of the initial data 
fields, although inadequate resolution in the numerical 
model could also explain the biases. 

Since the resources required to increase high quality 
data or to run sophisticated numerical models are expensive 
and time consuming, alternative methods for improving model 
forecasts should be considered. S ta ti s ti ca 1 - dynami ca 1 
schemes for predicting tropical storm motion (Renard et al, 
1972) rely heavily on current storm motion. If position 
reports of tropical storms are not timely, then the statisti- 
cal schemes may be impaired. An approach which circumvents 
many of these problems is to use the model forecast positions 



14 



themselves to statistically modify the predicted track posi- 
tions. It is the primary objective of this thesis to de- 
velop and evaluate statistical regression equations for 
adjusting dynamically predicted storm tracks from the Tropi- 
cal Cyclone Model (TCM) used by the U.S. Navy Fleet Numerical 
Weather Central (FNWC). 

The basic model used for these experiments is the primi- 
tive-equation, three-layer, tropical cyclone model developed 
by Elsberry and Harrison (1971) and Harrison (1973). Although 
this model is capable of triply-nested operation (Harrison, 
1973), results from Ley and Elsberry (1976) show that the 
coarse and nested grids produced nearly identical results in 
a selected case study based on hand-analyzed data. In 1975 
FNWC adapted a coarse-mesh (2°), three-layer, dry version of 
this model for tropical cyclone prediction. Preliminary 
results with operational data were presented by Hinsman 
(1977) for 1975 and 1976 data. For the 1977 typhoon season, 
the model was modified to include a biasing technique sug- 
gested by Shewchuk and Elsberry (1978). The use of the 
forecast stream function field to determine the cyclone 
forecast positions reduced errors in the relative vorticity 
tracking used prior to this modification (Shewchuk, 1977). 

The current operational version of the model has boundary 
conditions which are insulated, free-slip walls on the 
north and south and cyclic on the east and west. Since it 
was suspected that these boundary conditions could adversely 
affect the forecast storm track, Hodur and Burk (1978) in- 
corporated one-way interactive boundaries. Substantial 
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improvement in storm track prediction was made compared to 
the previous channel version of the model. The open 
boundary version of the TCM was used in the experiments 
described herein. 
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II. BACKWARD INTEGRATION 



A. THE MODEL 

A detailed description of the Tropical Cyclone Model 
(TCM) is presented in Appendix A. The TCM is a three-layer, 
p r i mi t i ve- equa t i o n model in pressure coordinates. It is 
an open boundary model with one-way interactive boundary 
conditions on the north and south and cyclic conditions on 
the east and west boundaries. 

B. MODEL MODIFICATIONS 

All forward integration of the model was carried out 
using the version of the TCM described by Hodur and Burk 
(1978). The backward integration was carried out using a 
negative time step of -600 seconds unless the northern 
boundary location was greater than 40°N, In the latter 
case, the negative time step was reduced to -450 seconds 
after the storm moved a sufficient distance to the north. 

In the forward integration, heat was added to the storm 
center as defined by a minimum wind at 1000 mb. The pur- 
pose of the heating function is to counteract the dispersion 
of the vortex due to the finite differencing (Ley and 
Elsber ry, 1976). If heat had been added during the back- 
ward integration, this would have contributed to a better 
definition of the storm center. However, this presents a 
physically unrealistic situation for a typhoon moving back- 
ward in time. For this reason, the heating function was 
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set equal to zero for all backward integrations, 
heating was rejected because it led to premature 
of the storm, which became difficult to track. 



Negative 
di spersion 
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Ill, APPROACH TO TRACK MODIFICATION 



Two types of regression equations were tested here. One 
type will have predictors based only on forward integration 
of the TCM, while the other set of predictors will be based 
on both forward and backward integration of the TCM. The 
open boundary version of the TCM (Hodur and Burk, 1977; Hodur 
and Burk, 1978) was run with 46 cases (9 storms) from 1975 
and 1976. All cases were based on o pe ra ti o na 1 1 y- a na 1 y ze d 
data obtained from FNWC. The regression equations developed 
from this sample were then tested against independent cases 
in 1977 and 1978. A second set of regression equations based 
only on forward integration of the TCM are derived using the 
combined data sets of 1975-78. 

The final set of regression equations are based on forward 
integration of the model to 72 hours, as well as backward 
integration for 36 hours. The backward integration should 
reveal the effects of systematic model and data errors. 
Fundamental assumptions are that the mode 1 - rel ated errors 
tend to be systematic, and adjustment for data errors is 
possible where, in the absence of observations, the initial 
analysis reverts to the east-west flow appropriate to the 
climatology. Backward integration will increase the number 
of predictors available to explain the variance between TCM 
forecasts and best track positions. This should lead to 
improved regression equations with a higher explained variance. 
The predictors using ba ckwa rd- i nteg ra te d positions can be 
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used to adjust the forward - 1 integrated positions using the 
same i ni ti al data . 

Both the forward- and backward- predi cted tracks were 
compared to the best tracks at their corresponding times. 
Next, regression equations to adjust the predicted track to 
the actual track using stepwise regression were generated 
with the Statistical Package for the Social Sciences (SPSS). 
The regression equations were developed to make corrections 
every 12 hours up to 72 hours. If a TCM run was incomplete, 
alternate sets of equations were derived according to the 
length of the TCM run. The options considered were 36-, 48- 
60-, and 72-hour TCM runs, with and without backward integra 
t i on . 

If successful, the advantage of this approach would be 
that use of simple regression equations would require much 
less computer time than more sophisticated dynamic models. 
Thus, it may be possible to produce tracks that are more 
accurate than -the ordinary open boundary TCM without a 
large increase in computer resources. An operational advan- 
tage of this approach is that no warning positions are 
required for computation of any predictors in the regression 
equations. The method and results of this approach to storm 
track prediction are presented in the following sections. 
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IV. DEVELOPMENT OF THE REGRESSION EQUATIONS 



A . THE PREDICTORS AND PREDICTANDS 

Using TCM forecast runs versus best tracks (Annual Typhoon 
Reports, 1975-77), 12 predictands were derived by computing 
the east-west and north-south differences between positions at 
corresponding times (see Fig. 1). Storm positions at each 
forecast time are adjusted by two regression equations, one 
for the east-west direction and the other for the north-south 
direction. Thus, for a 72-hour TCM forecast run, a total of 
12 regression equations would be used to modify storm tracks 
in 12-hour increments from 12 to 72 hours. 

Predictors used in these equations were model -predi cted 
displacement and velocity, broken into components along the 
east- we st and north-south directions. The Julian day, lati- 
tude and longitude of the initial position of each TCM fore- 
cast run are also included as predictors in the regression 
equations. A schematic illustration of the intervals over 
which the predictors were calculated for the forward integra- 
tion runs is shown in Fig. 2. A complete list of predictands 
and predictors, with the times for which they were computed, 
appears in Table I. 

B. METHOD OF EQUATION DERIVATION 

The regression equations were derived using the Statisti- 
cal Package for the Social Sciences (Nie, et al, 1975). 

Cases with missing values for predictors or predictands were 
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Figure 1. Depiction of the model errors (i.e., the predic- 
tands), which are the difference between the best 
track and forecast positions, which are shown 
above as AX^> AY^ AX^, ^ Y 24’ etc - 
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Figure 2 . Depiction of intervals over which displacements 
and average speeds were computed using only for- 
ward integration of the open boundary TCM. Un- 
modified TCM forecast positions (0) are shown at 
12 hour intervals. 



23 



TABLE I 



Predictors/predictands used to develop regression equations 
for typhoon track modification based only on forward inte- 
gration of the open boundary TCM. 

1. Predictands: AX, AY 

Times at which predictands are computed: 

12, 13, 26, 38, 60, 72 hrs 

2. Predictors: AX, AY, u, v 

Time intervals over which each predictor was computed: 
00-12 , 12-24 , 24- 36 , 36-48, 48-60 , 60-72, 

00-24, 12-36, 24-48, 36-60, 48-72, 00-36, 

00-48, 00-60, 00-72 hrs 

3. Initial Position Predictors: Julian Day, Latitude, 

Longi tude 

Times : 

Initialization time of a given TCM forecast run 
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automatically eliminated from all calculations. Such cases 
arise because the tracking routine is not always able to 
follow the storm center throughout the 72-hour interval. In 
the forward integration tests of Hodur and Burk (1978) with 
the open-boundary conditions, only 28 of the 46 cases extended 
to 72 hours. 

To avoid the problem of uncomputable predictors when 
the duration of TCM runs was less than 72 hours, it was 
decided to derive alternate sets of equations based on the 
duration of the TCM forecast run. The forecast lengths con- 
sidered were 36, 48, 60 and 72 hours. For example, if the 
TCM produced only a 48-hour storm track forecast, no regres- 
sion adjusted track would extend behond 48 hours, and only 
predictors in the range 00-48 hours were considered when 
deriving the regression equations to be applied to a 48-hour 
TCM run. 

All predictors listed in Table I were considered for in- 
clusion in each regression equation. Selection of predictors 
was stopped when the next variable in the stepwise regression 
explained less than 1 % of the variance. 

C. FIRST EQUATION SET BASED ON FORWARD INTEGRATION 

The 1975-76 cases, as well as those in 1977 and 1978 to 
be discussed later, are based on cases with relatively well- 
developed storms which are more readily modeled by the coarse- 
mesh TCM. It should be noted that the operational version of 
the TCM is used only for storms exceeding 50 knots. Inclusion 
of weaker tropical storms would likely increase the variance 
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between the forecast and best track positions, thus making 
it more difficult to derive stable regression equations. A 
sample equation for the regression adjustment along the x- 
axis at 7 2 hours is shown below: 

DXER72=-537. 0151+39. 2088(XXLAT)+68.1605(VX1224) 

- 50. 1251 (VX0024)-43.8154(VY0048) 

- 1 . 4 95 6 ( DX6 0 72 ) + 0 . 6 102 ( JULDAY ) 

Velocity was the most frequently -sel ected predictor. This 
result is not too surprising if one recalls that a common 
fault of dynamic models is prediction of motion which is too 
slow (Hovermale et al, 1976; Ley and Elsberry, 1976), The 
version of the TCM (Hodur and Burk, 1973) used in these ex- 
periments is no exception. Predictors with their associated 
time intervals beginning, ending or overlapping the valid 
time of a given regression equation were often selected. In 
general, this suggests a sort of "statistical extrapolation" 
in which past, present and future forecast motions are used 
to correct the forecast track. 

The average explained variance of the regression equations 
appears in Table II. This parameter measures the strength 
of the linear relationship between the multi-linear regres- 
sion equation value and the observed value. The amount of 
explained variance in Table II generally decreases as the 
TCM runs become shorter in duration and the number of avail- 
able predictors is thereby reduced. In general, longer TCM 
forecast runs permitted use of regression equations which 
corrected for more of the variance from the best track. 
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TABLE II 



Average explained variance of the regression equations. 



Forward Integration (1975-76 Cases) 



72 Hr TCM Run 81.6% 

60 Hr TCM Run 77.6% 

48 Hr TCM Run 75.4% 

36 Hr TCM Run 68.4% 

Forward Integration (1975-78 Cases)' 

72 Hr TCM Run 51.5% 

60 Hr TCM Run 46.5% 

48 Hr TCM Run 47.6% 

36 Hr TCM Run 43.2% 

Forward and Backward Integration (1977-78 Cases) 

72 Hr TCM Run 86.2% 

60 Hr TCM Run 84.3% 

48 Hr TCM Run 83.8% 

36 Hr TCM Run 82.3% 
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Shortcomings (see Fig, 3) of the one-way (OW) interac- 
tive boundary version of the TCM are that the forecast 
track typically runs to the left of the storm path and the 
velocities are generally too slow (Hodur and Kurk, 1978). 

These systematic errors of the TCM seem to contribute 
significantly to the large amount of explained variance in 
the regress ion-produced typhoon positions (see Table II). 

It should be recalled that the 1975 typhoon season, 
and to a lesser extent the 1976 season, experienced a high 
frequency of storm recurvature as well as a large number of 
storms which tracked northward. This bias was present in 
the data sample used to derive the regression equations and 
had considerable impact on the selection of predictors and 
the computation of constants and coefficients. 

D. SECOND EQUATION SET BASED ON FORWARD INTEGRATION 

To increase the number of typhoon cases used to derive 
the regression equations, TCM forecast runs from 1977 and 
1978 were included in the data set. These TCM forecasts 
were originally used to make an independent test of the 197 
76 regression equations. However, the majority of the 
typhoons in 1975 and many in 1976 had significant recurvature. 
The inclusion of the 1977-78 storm tracks should make the 
sample more representative of the various tracks found in 
the western North Pacific Ocean area. 

Because some TCM runs did not extend to 72 hours, only 
61 of 90 cases were available for derivation of the regres- 
sion equations. As with the first set of equations based 
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in 




Figure 3. Forecast tracks of Typhoon Marie produced by 
OW model compared to the 8-13 April 1976 bes 
track positions. Each circle represents a 1 
increment (after Hodur and Burk, 1973). 
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the 

-hour 



on the forward integration of the TCM, alternate sets of 
equations were derived based on 72-hour, 60-hour, 43-hour 
and 36-hour TCM runs. When less than \% of the variance was 
explained by a variable, selection of predictors was halted. 
The complete set of regression equations is listed in 
Appendix B. 

The expanded data sample (61 cases) contained a greater 
variety of storm tracks and reduced the recurvature bias of 
the 1975-76 typhoon cases. Storm tracks in the 1975-78 
data set used in these experiments are characterized by four 
general categories: westward, northwestward, northward, 
and recurving paths. With this greater variety of storm 
paths, the amount of explained variance of the regression 
equations (see Table II) dropped commens ura tel y . For example, 
the regression equations for a 72-hour TCM run incurred 
approximately a 30 percent decrease in explained variance. 
With the exception of a small variation at 48 hours in Table 
II, the general trend was again a reduction of explained 
variance as the duration of the TCM forecast runs decreased 
to 36 hours. As discussed previously, this characteristic 
is attributed to having fewer predictors available to explain 
the variance when the TCM forecast is of shorter duration. 
These results reaffirm that this method of adjusting TCM 
forecasts is at its best when the TCM forecast extends to 
72 hours, rather than a shorter forecast interval. 

Velocity was again the most frequently selected type of 
predictor, indicating an attempt to compensate for the over- 
all slowness of the TCM. However, there was a better 
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balance of velocity and displacement predictors in the 1 975- 
78 equations as compared to the 1975-76 equations. 

E. EQUATION SET BASED ON FORWARD AND BACKWARD INTEGRATIONS 

The third set of regression equations was based on storm 
positions derived from both backward and forward integration 
of the TCM. Additional velocity and displacement predictors 
based on backwards integration of the TCM as indicated in 
Fig. 4 and Table III were added to the data set. These pre- 
dictors are computed in the same manner as those based on 
forward integration. All predictors listed in Table I and 
Table III were considered for inclusion in each regression 
equation. Regression equations were again developed for TCM 
forecast runs of 36, 48, 60 and 72 hours duration. 

The regression equations which included ba ckwa rd- i n te- 
grated positions had the highest average explained variance 
(see Table II). As TCM forecast runs become shorter in 
duration, the amount of explained variance in the equations 
decreases due to fewer predictors being available (see 
Table II). The complete set of regression equations is 
listed in Appendix B. 

The regression equations took advantage of the systematic 
errors inherent in the initial fields and in the numerical 
model. The most favored predictors from the backward inte- 
gration occurred in the interval from -12 to 00 hours. For 
the 12 equations used when the TCM was integrated forward 
to 72 hours and backward to -36 hours, predictors in the 
interval -12 to 00 hours appeared in 9 equations. The 
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72 o 




Figure 4. Depiction of intervals over which displacements 
and average speeds were computed using both for- 
ward and backward integration of the open bound- 
ary TCM . Unmodified TCM forecast positions (0) 
are shown at 12-hour intervals. 
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TABLE III 



Additional predictors used to develop regression equations 
for typhoon track modification based on the backward inte- 
gration of the open boundary TCM. 



Predictors: AX, AY, u, v 

Times = 00-M12, M12-M24, M24-M36, 12-M12, 

00-M24, M12-M36, 00-M36, 24-M24 hrs 



where 00 = i n i ti al ti me 

M24 = minus 24 hours from initial time 
12 = plus 12 hours from initial time 
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velocity along the y-axis from -12 to 00 hours appeared In 7 
of the equations, and in each case it explained the most 
variance of any predictor in the equation. Approximately 
35 percent of all the predictors in these 12 regression 
equations were computed from ba ckwa rd- i n te gr a ted storm 
positions. Velocity was the most frequently selected pre- 
dictor from all forecast intervals. As will be indicated 
in later examples, the regression equations appeared to be 
compensating for the generally slow motion of the TCM. It 
should be noted here that these regression equations and 
the sample cases which follow were based on a limited data 
set of 31 storm cases from 1977 and 1978, 
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V. RESULTS 



A. FIRST EQUATION SET BASED ON FORWARD INTEGRATION 

The primary purpose of the tests with the 1975-76 cases 
was to make a preliminary evaluation of the regression 
equations which were derived using this data set. Because 
this data set was used to derive the regression equations, 
these experiments are referred to as a dependent test. 

A sample of adjusted storm positions based on the first 
equation set is shown in Figs. 5-7. Results with Typhoon 
June shown in Fig, 5 and Typhoon Marie depicted in Fig. 6 
were very encouraging, Typhoon Marie is the same storm 
which appears in Fig. 4, Storm recurvature was accurately 
depicted, and storm velocity was markedly improved in these 
two storms. Note, however, that the velocities in the 
regression tracks are generally greater than the best track 
velocities. This characteristic of the regression equations 
appears to be an attempt to compensate for the slowness of 
the TCM, In the majority of cases in this sample, the 
storms were excessively accelerated due to this feature of 
the equations. Since the regression equations were based 
primarily on predictors derived from TCM forecasts, it is 
clear that the goodness of the regression-adjusted storm 
track is dependent on the quality of the TCM forecast itself. 
The adjusted track for Typhoon Rita (August 1975) shown in 
Fig. 7 illustrates how a poor TCM forecast will lead to an 
extremely radical regression adjustment. This example 
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BEST TRACK 
T CM FORECAST 
REGRESSION ADJUSTMENT 
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Figure 5. Unmodified TCM forecast tracks and regression- 
adjusted tracks of Typhoon June (18 Nov 75, 00 
GMT) compared to best track positions. The in 
terval between adjacent positions is 12 hours. 
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Figure 6. Typhoon Marie (10 Apr 76, 00 GMT), otherwise the 
same as Fig. 5. 



45. 



2: 



CO 

3 



Q 
I— <X 
CO 




» - 1 ! : * - - f - i i I : I 



o 

=r 



lo o 

on cn 



a 

o 



o 

LO 

cn 




LU 



CD 



CD 

o 



(N) sanmtn 



38 



Figure 7. Typhoon Rita (21 Aug 75, 12 GMT), otherwise the 
same as Fig. 5. 



indicates how small TCM velocities, as well as radical 
changes in the forecast storm track, will result in a 
regression-adjusted track with extremely high velocities 
and unrealistic variations in the storm path. This behavior 
of the regression equations seems to indicate that obviously 
erroneous regression positions can be used as a basis to 
reject TCM forecast positions as well. 

The mean forecast errors of the 1975-76 dependent test 
sample are listed in Table IV. The sample improved on the 
forecast errors of both JTWC and the unmodified TCM fore- 
casts for 1975-76. At 72 hours with a sample size of 23, the 
regression equation errors were also less than the U.S. 

Navy 7th Fleet goal of 150 nm. However, it is not expected 
that this pattern would be repeated in an independent sample 
of typhoon cases . 

B. INDEPENDENT TEST OF FIRST EQUATION SET 

An independent sample of 44 cases from the 1977-73 typhoon 
seasons was then used to test the regression equations de- 
rived from the 1975-76 cases discussed in the previous 
section. As indicated by the statistics in Table V, the 
results of the independent test were generally very poor. 

In each category the modified TCM tracks were worse than the 
unmodified tracks. The large errors are attributed to the 
unstable nature of the regression equations from the small 
and very homogeneous sample of anomalous storm tracks used to 
derive the regression equations. This led to excessive 
velocities and erratic tracks in the regression-adjusted 
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Mean forecast errors (nm) based on the best tracks for 1975 and 
1976. Errors incurred by the unmodified TCM and the regression- 
modified TCM are shown with the combined 1975-76 JTWC errors and 
the U.S. Navy 7th Fleet goal for forecast errors. 
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Mean forecast errors (nm) based on the best tracks for 1977 and 
1978. Errors incurred by the unmodified TCM and the regression- 
modified TCM are shown with dependent test results for 1975-76 
and the U.S. Navy 7th -Fleet goal for forecast errors. The same 
regression equations were used in the modified TCM positions for 
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storm positions for 1977-78, When the TCM forecasts had 
large deviations from the best track, the regression positions 
incurred commensu ra te 1 y larger errors. It was clear that 
the number of TCM forecast cases had to be increased as much 
as possible to increase the stability of the regression 
equations. 

The regression equations used in these tests did not 
improve storm positions by taking a-dvantage of the systematic 
errors of the TCM, which is often too slow and to the left 
of the best track. Typhoon Lola is depicted in Fig. 8 and 
indicates slight corrective shifting of the regression posi- 
tions to the right of the TCM forecast track. However, over- 
compensation is again evident in the velocity components. 

More typical of the independent test results is the adjusted 
track for Typhoon Lucy shown in Fig. 9. Velocity errors were 
large and the track was erratic. In many cases in the inde- 
pendent test, an erratic track, such as the regression- 
adjusted path -for Typhoon Lucy, occurred when the unmodified 
TCM forecast track also incurred large errors. The single 
redeeming feature of these 1975-76 equations may be that 
highly erratic regression positions could be indicative of a 
poor TCM forecast. This information may be useful in leading 
the typhoon forecaster to reject both the modified and the 
unmodified TCM guidance in making his decision. 

C. SECOND EQUATION SET BASED ON FORWARD INTEGRATION 

The second set of regression equations was obtained from 
the combined storm cases in the 1975-76 and 1977-78 samples. 
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Figure 8. Typhoon Lola (29 Sep 78, 00 GMT), otherwise the 
same as Fig. 5. 
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Figure 9. Typhoon Lucy (03 Dec 77, 00 GMT), otherwise the 
same as Fig. 5. 



Because this exhausts the sample, these equations can only 
be evaluated as a dependent test. 

This second set of equations, with a larger sample of 
61, was more stable overall than the first set of regression 
equations. While the regression-adjusted positions still 
increase the velocity excessively, the apparent result is 
that this overcompensation is decreased in the 1975-78 re- 
gression equations versus those for' 1975-76, 

The five typhoon cases shown previously in Fig. 5 through 
Fig. 9 appear again in Fig. 10 through Fig, 14, but this time 
illustrating the regression-adjusted positions of the 1975- 
78 equations. The track for Typhoon June shown in Fig. 10 
was not adjusted as well in this test, but it serves to 
illustrate that the regression adjustments cannot always com- 
pensate for the tendency of the TCM to not predict recurvature. 
Presumably if the TCM track had been more northwesterly, the 
adjustment would have been toward more recurvature. The 
track for Typhoon Marie depicted in Fig. 11 still indicates 
a reasonable adjustment for recurvature, but has incurred a 
marked decrease in velocity adjustment. Although the adjusted 
speed of movement for Typhoon Rita shown in Fig. 12 is slower 
than that on Fig. 7, the track is not significantly changed 
from the previous result. However, this case again suggests 
that a radical and obviously erroneous regression adjustment 
may serve as a basis to reject the TCM forecast also. Typhoon 
Lola (September 1978), illustrated in Fig, 13, has a much 
improved regression adjustment of the storm track compared 
with Fig. 8. The velocity corrections in this case are much 
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Figure 10. Typhoon June (18 Nov 75, 00 GMT), otherwise the 
same as Fig. 5. 
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Figure 11. Typhoon Marie (10 Apr 76, 00 GMT), otherwise the 
same as Fig. 5. 
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Figure 12. Typhoon Rita (21 Aug 75, 12 GMT), otherwise the 
same as Fig. 5. 
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Figure 13. Typhoon Lola (29 Sep 78, 00 GMT), otherwise 
same as Fig. 5. 
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Figure 14. Typhoon Lucy (03 Dec 77, 00 GMT), otherwise the 
same as Fi.g. 5. 



closer to the actual velocities which appear in the best 
track positions. The adjusted track for Typhoon Lucy is 
shown in Fig. 11 and is to be compared with Fig. 9. The 
regression-adjusted positions are still poor, but not 
nearly as erratic. Note that the regression-modified posi- 
tions do suggest some recurvature, although i nacc ura tel y , 
and that the decreased velocities are much more realistic 
with the 1975-78 regression equations. 

Generally, the results of these experiments suggest 
that a further increase in the number of TCM forecasts in 
the data sample would lead to a further improvement in the 
regression-modified storm positions. Overall, the errors 
incurred by the regression-adjusted tracks were less than 
those of the unmodified TCM storm path (see Table VI). The 
improvements are especially noteworthy at 48 hours and 72 
hours . 

In forecasting typhoons which tracked westward, the un- 
modified TCM forecast positions usually proved to be very 
accurate, with the regression-modified path seldom improving 
on the TCM. One case (Typhoon Rita, October 1978) in which 
the regression equations improved the forecast is shown in 
Fig. 15. The improvement is the result of an increase in 
storm velocity by the regression equations. Several other 
examples of enhanced track predictions are shown in Fig. 16 
through Fig. 18. These examples show some improvement in 
direction, but the improvement in speed of movement is 
especially rewarding. 
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Mean forecast errors (nin) based on the best tracks for 1975 
through 1978. Errors incurred by the unmodified TCM and the 
regression-modified TCM are shown with the U.S. Navy 7th 
Fleet goal for forecast errors. 
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Figure 15. Typhoon Rita (23 Oct 78, 00 GMT), otherwise the 
same as Fig. 5 . 
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Figure 16. Typhoon Phyllis (14 Aug 75, 12 GMT), otherwise 
the same as Fig. 5. 
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Figure 17. Typhoon Babe (06 Sep 77, 00 GMT), otherwise the 
same as Fig. 5. 
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Figure 18. Typhoon Pamela (22 May 76, 00 GMT), otherwise 
tite same as Fig. 5 . 



0. EQUATION SET BASED ON FORWARD AND BACKWARD INTEGRATION 

These experiments with the 1977-73 typhoon cases were 
performed to test the regression equations derived from the 
1977-78 data set which included backward-integrated posi- 
tions to -36 hours. As the complete set of backward tracks 
was used in the derivation of the regression equations, 
these experiments refer only to the dependent sample. 

The results of the dependent test with forward- and 
backward-integrated positions were the best of all experi- 
ments conducted during this research. This equation set in- 
curred the lowest mean errors overall (see Fig. 19). The 
regression-modified storm position errors are very close 
to meeting the U.S. Navy 7th Fleet error goal (see Table VII). 

Tracks for Typhoons Lola, Rita, and Babe shown previously 
in Fig. 13, 15, and 17, respectively, appear again in Fig. 20 
through Fig. 22 with adjusted positions based on forward- 
backward integration. Typhoon Lola (September 1978), shown 
in Fig. 20, is’ more erratic than the same case (Fig. 13) 
using regression equations based only on forward integration. 
However, the regression positions that are based on forward- 
backward integration are fluctuating on either side of the 
best track. The variation about the best track is attributed 
primarily to the small size of the data set (31 cases) used 
to derive the equations. This "saw-tooth" variation about 
the best track is characteristic of the behavior storms which 
tracked westward in this experiment. As in previous tests, 
it proved difficult to improve on the unmodified TCM forecast 
positions when typhoons were tra eking approximately westward. 
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FOREC ST TIME(HRS) 



Figure 19. Mean track error (nm) for the 1975-76 dependent 
test (0), the 1977-78 independent test (a), the 
1975-78 dependent test (•), and the 1977-78 
dependent test C a ) • The 1977-78 dependent test 
(a) was the only experiment which included back- 
ward integration. 
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Figure 20. Typhoon Lola (29 Sep 78, 00 GMT), otherwise the 
same as Fig. 5 . 
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Figure 22. Typhoon Babe (06 Sep 77, 00 GMT), otherwise the 
same as Fig. 5 . 



The adjusted track of Typhoon Rita in Fig. 21 exhibits a 
similar behavior in another westward tracking storm, but 
does very well in the longitudinal positioning of the storm. 
Adjustment of the track of Typhoon Babe, shown in Fig. 22, 
has improved the velocity forecast as well as has indicated 
some recurvature which the TCM missed (see also Fig. 17). 

Typhoon Phyllis is illustrated in Fig, 23. It provides 
an excellent forecast of recurvature not achieved by the 
unmodified TCM forecast or by the regression equations based 
only on forward integration of the TCM. The velocities in- 
dicate only a slight tendency to overcompensate for the 
slower TCM forecast. Recurvature of Typhoon Gi 1 da , depicted 
in Fig. 24, is correctly predicted, but the velocities are 
excessive. Adjustment of the track of Typhoon Wendy as in 
Fig. 25 showed slight improvement over the unmodified TCM 
positions. The interesting feature in the regression-adjusted 
path is the indication that the storm would reverse its 
direction. Within 12 hours of the indicated time, Typhoon 
Wendy did change direction in this manner, but with somewhat 
smaller displacements. 
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Figure 23. Typhoon Phyllis (19 Oct 78, 00 GMT), otherwise 
the same as Fig. 5. 
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Figure 24. Typhoon Giltla (05 Oct 77, 12 GMT), otherwise the 
same as Fig. 5. 
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Figure 25. Typhoon Wendy (27 Jul 78, 12 GMT), otherwise the 
same as Fig. 5. 



VI. CONCLUSIONS 



A three-layer, primitive-equation model (Hodur and Burk, 
1978) with one-way interactive boundaries is being tested 
at the Naval Environmental Prediction Research Facility 
(NEPRF). The objective of this research was to generate 
statistical regression equations to adjust the TCM- p red i c ted 
tracks towards the best tracks. This approach is based on 
the assumption that it is possible to adjust for systematic 
model and data errors. Development consisted of deriving 
three sets of regression equations, with two sets based only 
on forward integration of the TCM, while the last set con- 
tained predictors based on both forward and backward inte- 
gration. All TCM runs were based on operationally analyzed 
data from FNWC. Ti me- de p en den t boundary conditions provided 
to the TCM were derived from analyzed (rather than predicted) 
fields, as was the case in Hodur and Burk (1978). 

The most notable improvements occured with the regression 
equations containing predictors based on backward integration 
of the TCM. The equations with predictors having both for- 
ward- and backwa rd- i ntegrated positions explained the greatest 
amount of variance of any set of regression equations. Ad- 
justments to the TCM tracks at 12-hour intervals resulted 
in predictions in the 1977-78 dependent test that nearly 
met the 7th Fleet forecast error goals. The selection of the 
velocity along the y-axis in the interval -12 to 00 hours as 
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the predictor explaining the most variance suggests strongly 
that systematic errors at or near initialization time were 
used advantageously to adjust TCM track forecasts. It is 
noteworthy to observe that TCM forecasts with all data samples 
used in these experiments incurred systematic errors (see 
Fig. 26). 

It should be noted that the weak link in the regression 
adjustment scheme is the small sample size for derivation of 
the regression equations. With the exception of the 1975-76 
equations, the regression coefficients have not been tested 
with an independent sample of storms. Application of the 
statistical equations in the future may not produce comparable 
results due to the relatively small number of cases used to 
derive the equations. Likewise some years have persistently 
anomalous storm tracks, as in the first sample (1975-76) 
treated here. It is expected that some increase in the 
stability of the regression equations could be achieved by 
increasing the'number of TCM forecast cases used to derive 
the equations. This is especially true for the backward 
integration set. 

In all of the experiments, the unmodified TCM forecasts 
were difficult to beat in the case of westward propagating 
storms. However, occasional improvement in velocity predic- 
tion was noted. The typhoons used in these tests had storm 
tracks in several general categories, such as westward, south- 
east- to- no rthwe s twar d , northward, or recurving. Since the 
regression equations were derived from storm cases involving 
various combinations of these track types, it was felt that 



AY ERROR (N Ml) 



AX ERROR (N Ml) 




Figure 26. Mean errors of the TCM relative to the best track 
positions indicate systematic error. The AX 
versus the AY error is shown for the 1975-76 cases 
(a), the 1977-78 cases (a), and the combined cases 
for 1975-78 (0) . 
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these widely varying tracks and subsequent errors contri- 
buted to failures of the regression equations to make proper 
track adjustments in some cases. A possible solution might 
be to vastly increase the TCM case sample size and then 
subdivide the cases into "regression analogs" according to 
the direction of typhoon propagation. The next step would 
be to derive a set of regression equations based on each 
subset of storm tracks. It could be expected that each sub- 
set of equations would have inherent adjustment biases 
towards the type of track from which they were derived. Any 
or all of these subsets of equations could be applied to each 
unmodified TCM forecast. 

On the basis of the sample examined here, it appears that 
adjustment for recurvature was handled most effectively by 
the regression equations which included backward- integrated 
tracks. It is suggested that regression equations based on 
an expanded sample of TCM cases including backward integration 
should be derived to produce better adjustments and statis- 
tical equations with greater stability. 
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APPENDIX A: THE TROPICAL CYCLONE MODEL 



A. THE FORECAST MODEL 

The coarse-mesh version of the primitive equation model 
developed by Elsber ry and Harrison (1971), Harrison (1973), 
Ley and Elsberry (1976), and Hodur and Burk (1978) was used 
as the basis for these experiments. The model is a coarse 
mesh (2°), three-layer, dry model with one-way interactive 
boundaries on the north and south and cyclic east-west 
boundaries. The equations used in the model are: 
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L(S) - H^tf) ♦ |y(^f)] - Ip <«s) 
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L(S) represen ts the flux divergence of any scalar quant i ty S. 
The meteorological symbols used above can be found in the 
"List of Symbols". 

A sufficient condition for the linear computational sta- 
bility of the solution for two-dimensional equations govern- 
ing simple wave motion is: 

^<•707. (A-7) 

where 

C = phase velocity of the fastest gravity wave 
At = time increment 
Ax = horizontal grid increment 
Computational stability requires a maximum time step of 450 
seconds for this model. Ley (1975) achieved an increase to 
800 seconds by time averaging the pressure gradient term of 
the momentum equations. Shewchuk (1977) used a 600 second 
time step for testing his 1975 cases as the relocatable grid 
was extended northward where Ax is reduced. In the TCM ver- 
sion (Hodur and Burk, 1978) used for these experiments, the 
time step was 600 seconds unless the northern boundary exceeded 
40°N, in which case the time step was 450 seconds. 

The initial step was forward in time, while the leapfrog 
time differencing scheme was used in all subsequent iterations. 
Friction was neglected, and the consequent storm motion was 
primarily the result of advective processes and heating. A 
Bessel interpolator was used to locate the minimum wind at 
1000 MB. Then latent heating was simulated by adding heat to 
a horizontal 7 x 7 grid centered on the minimum wind, which 
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was not necessarily on a grid point. Weighting of the heat- 
ing function in the vertical was 0.3, 1.0, and 0.3 at 350, 

550, and 250 MB respectively. In the horizontal, the effect 
of heating is smoothed out by a Cress man weighting function 
(Haltiner, 1971) and results in a less erratic storm track. 

B. THE GRID 

The forecasts are carried out on a uniform, coarse-mesh 

/ 

(2°) Mercator grid true at 22.5°N. The horizontal grid inter- 
val was 205.8 km. The domain consisted of 32 points east- 
west and 24 points north-south. The grid was oriented so 
that each storm was initially located southeast of the center 
of the grid. The vertical distribution of variables is shown 
in Fig. A - 1 . Although the variables are staggered in the 
vertical, staggering on the horizontal grid was not used in 
these experiments. The 850 mb winds were used to compute the 
advective terms at 1000 mb. 

C. BOUNDARY CONDITIONS 

Boundary conditions on the walls were after Hodur and 
Burk (1978). In the one-way interactive (0W) model, wind com- 
ponents normal to the boundaries are adjusted so there is no 
net divergence from the forecast domain. By integrating 

If = V n (A-8) 

around the boundary, the corresponding values are obtained 
(Hawkins and Rosenthal, 1965). Distance along the boundary 
is represented by s, while V n is the velocity component normal 
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Fig. A-l. Vertical distribution of dependent variables 
and pressure levels for the three-dimensiona 
model (after Harrison, 1973). 
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to the boundary with the positive direction being inward. 

Date from outside the forecast grid must continuously 
be incorporated into the boundaries. This is accomplished 
as described by Perkey and Kreitzberg (1976), and following 
their notation, the prediction equation for a variable X 
becomes 






j ) 



X p (i ,j)+W(i ,j) 



3X. 



m 



at 



A t + [ 1 - W ( i , j ) ] 

» j 



ax 

Tt 




At 



. j 



( A- 9 ) 



where the subscripts n and p indicate the new and the previous 
value of X. The subscript m refers to the model tendency, 
while Is is the large-scale tendency. The weighting function 
used in this model was somewhat different from that given by 
Perkey and Kreitzberg (1976). 

The following weighting function W(i,j) was constructed 
so that a minimum amount of noise was produced near the 
boundaries 



0.0 on the boundaries 
0.05 one grid row in from boundaries 
0.25 two grid rows in from boundaries 
0.45 three grid rows in from boundaries 
0.65 four grid rows in from boundaries 
0.85 five grid rows in from boundaries 
1.00 on all other interior points 
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The forecast fields produced must be filtered near the 
boundaries. The filter used in this model is 

FT = (l-a)F j + (f)(F j + 1 + F j _ 1 ) (A- 10) 

where F. is the filtered data at point j, F. is the unfiltered 
J J 

data at point j, and a represents the smoothing parameter. 

The filter was applied every 40 minutes to the six rows and 
columns nearest the boundaries with a = 0.5. This value of 
a yields a response function of 0 for a 2Ax wave. This filter 
is also applied over the entire grid to all the prognostic 
variables every three hours. In this manner the forecasts 
include large-scale tendencies from outside the forecast 
domain. 

The vertical velocity at the upper boundary is equal to 
zero and is calculated at levels 5, 3, and 1, Fig. A-l, 
through downward integration of Eq. (A-5). 

In an operational mode, the time dependent boundary condi- 
tions must be specified by a global forecast model. In these 
experiments, forecast fields were not available. Therefore, 
using a "perfect-prog" approach, analyzed fields taken every 
12 hours were used to specify the boundary values. 

D. INITIALIZATION 

Analyzed fields obtained from FNWC are used as input data 
for the model. Initialization of the model is accomplished 
by calculating non-di vergent winds from the stream function j. 
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The relative vorticity, the forcing function for the stream 
function, was obtained from the observed u and v components 







( A- 1 1 ) 



The normal component of the wind was adjusted so there is no 
net inflow or outflow in the domain. The stream function was 
found from the expression 




( A- 12 ) 



A direct solver (see Faulkner and Rosmond, 1976) rather than 
successive over-relaxation is used to efficiently solve the 
Poisson equations in the initialization process. The non- 
divergent wind components were then calculated through 



= -M 
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= M i! 

5 x 



(A-l 3) 



An appropriate balance between the mass and motion fields 
is achieved through partial differentiation of Eq. (A-l) with 
respect to x and Eq. ( A - 2 ) with respect to y. Addition of 
these equations and assuming that the time rate of change of 
divergence and gradients of the map scale factor can be neg- 
lected leads to a Poisson equation 



7 2 <t, = - J r (f 7 n.(u 4 ,)n- 



fy [L ( V ] + f ^ • 7T 1 '% T7 5 



( A- 1 4 ) 



which can be solved using direct methods. 
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APPENDIX B. THE REGRESSION EQUATIONS 



A. FORWARD INTEGRATION (1975-78 DATA) 



The regression equations which follow were produced by 
SPSS using only predictors derived from forward integration 
of the TCM. The data sample consisted of 90 cases, or runs 
of the TCM. 

1 . TCM Integrated to 72 HR 



VARIANCE 

REGRESSION EQUATIONS EXPLAINED 

DXER12=-110. 7609+ 1.5855 ( 0X1224 )+ 0 . 9 1 54 ( D Y2 4 36 ) 52.4% 

- 1 3 . 91 66 ( V X0024 ) + 4 . 7 87 5 ( X XLAT ) 

- 0. 3990 ( DY0012 )- 0 . 2 31 9 ( J UL DA Y ) 

+ 2 .1336 ( VX6072 ) 

DYER1 2= 33. 3375 + 18. 1 0 1 0 ( V Y 1 2 2 4 ) - 0 . 56 74 ( D Y00 36 ) 60.8% 

+ 4 . 4 54 7 ( XXL A T ) - 1 . 10 54 ( DX 3648 ) 

+ 17.0027 (VX 3660)- 3 . 4 1 88 ( V X 5 0 7 2 ) 

- 0. 3610(DY6072 )- 0 . 2 3 3 3 ( JU L DA Y ) 

DXER24 = -261 . 2548+ 0.0542 ( 0X2435 )+ 2 . 220 7 ( D Y24 36 ) 51 . 5% 

+ 11 . 2215(XXLAT) - 0 . 2442 ( DY4860 ) 

+ 13.0592(VY1224)-19.3594(VX0012) 

- 20.6092(VY0036)- 0 . 4907 ( J'JLDA Y ) 

+ 245. 7376(VX 0036)+ 6.1574(VX5 072 ) 

-289. 1865 (VX0048) +67. 2472 (VX3648) 

D Y E R 2 4 = 371.6372+ 4 . 7059 ( XXLAT ) -63.0716(VX6072) 54.9% 

+ 6. 3397(VY4872)-10.9032(VY0012) 

- 0.0693(DX1224)+18.3609(VX2436) 

- 1 .4547(XXL0N) - 4.9576(0X0048) 

+ 326.2109(VX0072)+ 1 . 51 70 ( DY0024 ) 

- 41 .8338 (VYOO 72) -4 9. 6258(7X4860) 

DXER36=-344. 7210+ 2 . 2 54 8 ( D Y 2 4 36 ) + 0.7540(0X1224) 45.0% 

+ 15. 3565(XXLAT) +20 . 9 332 ( V Y 1224 ) 

+ 0 . 8946 ( DX6072 ) - 1.0326(0X3643) 

- 0. 5088( JULDAY)-27.7637(VY0036) 

- 24. 1 1 94 ( V X00 1 2 ) + 2 7. 40 70 ( VX0036 ) 

DYER36= 260.5690- 0.9549(0X0024)- 2 . 9 2 5 5 ( D Y 2 4 36 ) 55.2% 

- 1 . 3287 (DX3648)- 0 . 5 7 72 ( D Y6 0 72 ) 

- 1 5 . 52 79 (VX6072 )+ 0 . 5636 ( DX0072 ) 

+ 0. 5077(DY0036) 
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4 3.3% 



DXER43=- 37.0420- 3 . 6 9 1 5 ( V Y2 44 8 ) + 2.9949(0X1236) 

- 96.8503(VX0048)+16. 3630(XXLAT) 

+ 1 . 8046 ( D Y 1 2 36 ) + 3.7131 (VX6072) 

+ 1 . 9662 ( DX2436 )- 2 . 4 9 88 ( X X L ON ) 

D YER48= 421.5829- 6 . 0 30 9 ( V XOO 1 2 ) - 2 . 2 2 9 1 ( D Y2 4 36 ) 55 . 5% 

- 20. 5254(VX0072 )- 0 . 50 54 ( DY60 72 ) 

DXER50= 97.7600- 9 . 30 39 ( V Y244 8 ) + 16 . 2 2 2 7 ( VX60 72 ) 37 . 1% 

+ 21 . 3398(XXLAT ) + 6.1324 ( 0X24 36 ) 

- 59.0393(VX2448)+ 1 . 9 7 5 3 ( D Y 1 2 36 ) 

- 22.6283(VX0012)- 3 . 7 7 82 ( XX LON ) 

+ 0. 72 77 ( DY60 72 ) 



0 YER50= 426. 2 759- 14. 02 0 7 ( V XOO 1 2 ) - 3.. 9702 (DY2436) 61.5% 

- 0. 132 7(0X2 44 8)- 16. 1361 ( V Y 4 8 7 2 ) 

+ 0.9047 (DY0036 )- 8 . 6 6 88 ( V X 6 0 72 ) 

OX ER 7 2 = 332.5937+ 7 . 1 4 91 ( D Y 24 36 ) + 20 . 2 1 35 ( V X60 72 ) 33.3% 

+ 33. 3952 (XXLAT) - 5 . 1 61 6 ( XX LON ) 

- 9 . 6 322 ( VYOO 12 ) - 1 . 0 394 ( J IJL DA Y ) 

- 3 1 . 326 2 ( V Y4860 ) + 3.7294 ( 0X2436 ) 

- 4 1 . 62 53 ( VX0048 ) 

0 Y ER 7 2 = 51 3. 8513- 19. 1 324 ( VXOO 1 2 ) - 34 . 76 30 ( V Y4872 ) 67 . 5% 

- 1 1 . 1924 (VX6072)- 3 . 0878 ( DY24 36 ) 

+ 32 . 4353 ( VYQ924 ) 



2 . TCM Integrated to 60 HR 
REGRESSION EQUATIONS 



,0 

VARIANCE 

EXPLAINED 



DXER12 = - 107.0123+ 1.6448( 0X1 224 )+ 1 . 1 26 5 ( D Y 2 4 36 ) 52.0% 

- 1 1 . 9798( VX0024)+ 4 . 3382 ( XXLAT ) 

- 0. 3577(0 Y 0012)- 0 . 2 1 5 2 ( J U LDA Y ) 

- 0.2897(DY4360) 



0 YE R 1 2 = 1 8. 042 5 + 15. 8288 ( V Y 1 2 24 ) - 0 . 5743 ( DY0036 ) 52.8% 

+ 3. 5755 ( XXLAT ) - 0.9042(0X3643) 

+ 8.4002(VX3560) 



DXER24 = - 294 . 5950+ 0.5243( 0X2436 )+ 2 . 380 3 ( 0 Y 2 4 36 ) 47.7% 

+ 1 0 . 6428 ( XXLAT ) - 0 . 6 9 73 ( DY4860 ) 

+ 1 5 . 351 7 (VY1224 )- 18 .4575(VX001 2) 

- 16.4373(VY0036)- 0 . 3287 ( JULOAY ) 

+ 24 . 486 5 ( VX0036 ) 

DYER24= 315.0132+ 3 . 8023 ( XXLAT ) - 1 . 2 2 7 1 ( D X 36 4 8 ) 37.9% 

- 1 . 0044 ( DY2436 ) + 7 . 6 2 86 ( V X 36 6 0 ) 

- 0 . 9660 ( XXLON ) 
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DXER36=-401 . 1934 + 2 . 34 6 1 ( D Y 24 36 ) + 2.7307 ( 0X1224 ) *1.1% 

+ 14 . 1 81 2 ( XXLAT ) +2 5 . 5 7 85 ( V Y 1 224 ) 

- 2 9 . 8364 ( V Y9060 ) + 2.0307(0X2430) 

- 70. 5401 (VX0048)+ 0.9644(0X2443) 

- 0. 32 35 ( JULDAY) 

0 Y E R 35 = 323.7302- 0.3010( 0X0024 )- 2 . 2 99 7 ( D Y24 36 ) 46.5% 

- 0.9573(0X3648) 

D YER48= 11. 4836- 10. 2403 ( V Y2448 ) + 2.9591 ( 0X 1236 ) 42.3% 

-127.1613(VX0048)+15.9606(XXLAT) 

+ 2 . 06 26 ( DY1 2 36 ) + 0.7195 ( 0X0060) 

+ 1.8145 ( 0X2436 )- 2 . 7430 ( XXLON ) 

D YER43= 41 7.1 723- 6 . 599 7 ( V XOO 1 2 ) - 2 . 44 56 ( D Y 2 4 36 ) 52 . 1% 

- 20 . 9459 ( V X0060 ) 

DXER60= 121 .4777 + 19. 3487( VY244S)+6 5. 1 2 9 9 ( VX2436 ) 33.3% 

+ 20 . 7 381 ( XXLAT ) + 2 . 2 9 1 0 ( D Y 1 2 36 ) 

- 23. 5426 (VX0012)-30. 3043(VX2443) 

- 4 . 0 5 3 7 ( XX LON ) - 2 4 . 33 7 3 ( V Y 366 0 ) 

0 YER60= 436. 1107- 12. 6669( VX0012) - 3 . 8492 ( DY2436 ) 59.1% 

- 0 . 6 3 1 6 ( DX2448)+23. 5083 ( VY0036 ) 

- 0 . 7 763 ( DY 4360 ) 



3 . TCil Integrated to 43 HR 
REGRESSION EQUATIONS 



<1 

/O 



VARIANCE 

EXPLAINED 



0 X ER 1 2 = -104.2041+ 1.7146 (0X1224 )+ 0 . 8 7 20 ( 0 Y 2 4 36 ) 50.7% 

- 14. 3293(V XOO 24 ) + 3 . 6437 ( XXLAT) 

- 0.3320(DY0012)- 0 . 1 56 2 ( J (JL DA Y ) 



DY ER1 2= 22. 7200 + 13. 7975 ( VY1224 ) - 0 . 5475 ( DYOO 36 ) 59.6 % 

+ 3 . 87 73 ( XX LAT ) - 0.0330( 0X3643 ) 

+ 15 . 9946 ( VX2436 )-21 . 5329( VX0043) 

- 0 . 3737 ( DY3648) 

DXER24 = -313.4117+ 1.5624( 0X2436 )+ 1 . 0949 ( 0 Y24 36 ) 46.5% 

+ 7 . 22 52 ( XXLAT ) +1 1 . 46 32 ( DX 1 2 24 ) 

-513. 1382 (VX0048)+ 0 . 7 74 6 ( 0 Y 1 2 24 ) 

+ 1 0 . 55 39 ( DX244S) + 9 . 850 2 ( DXOO 1 2 ) 



DYER24= 332.6320+ 4 . 7761 ( XXLAT ) - 0.8977(0X3648) 45.4% 

- 1 .6715(0Y2436)- 0 . 5969 ( DY0048) 

+ 29 .4431 ( VY1236 ) - 1 . 0 56 1 ( XXLON ) 

- 5.7302(VX0012)+15.0708(VX2436) 

+ 0.4376( D Y 00 1 2 ) - 8 . 6 5 7 5 ( V X 1 2 36 ) 
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DXER36 = -44 9 . 5619+ 1 . 2046 ( D Y 24 36 ) + 1 4 . 0084 ( DX 1 2 2 4 ) 
+ 1 1 . 1369 ( XXL AT ) + 16 . 31 31 ( VY 1224 ) 

+ 1 . 9 548 ( 0X2436 ) -6 13. 931 7 ( VX0048) 

+ 12.3598(0X2448)+ 11.3264(0X0012) 

DYER36= 323.7302- 0.3010(0X0024)- 2 . 2 9 9 7 ( D Y 2 4 36 ) 

- 0.9573(DX3548) 

DXER48=-512. 1396-1 3. 7846 (VY2 448)+ 3.2353(DX1236) 

- 91. 21 84 (VX 0048)+ 1 3 . 7796 ( XXLAT) 

+ 1 . 849 3 ( DY 1 2 36 ) + 1 . 7 2 1 7 ( D X 24 36 ) 

D Y ER4 8= 415.0658- 6 . 31 1 2 ( VX00 1 2 ) - 1 . 86 54 ( D Y 24 36 ) 

- 0 . 7522 ( DX2448 ) - 0 . 44 7 4 ( DY 2448 ) 

4 . TCM Integrated to 36 HR 
REGRESSION EQUATIONS 

DXER12=-104. 2041+ 1 . 7 1 46 ( DX 1 224 ) + 0 . 8720 ( DY 2436 ) 

- 14. 3293(VX0024)+ 3 . 6487 ( XXLAT ) 

- 0.3320 (DY0012 )- 0 . 1 562 ( JULDAY ) 

DYER12=- 5 .4141+18. 6 1 96 ( V Y 1 2 2 4 ) - 0 . 5 78 1 ( D Y 00 36 ) 
+ 4.4260(XXLAT) - 0.2211(0X1224) 

+ 12.1161 (VX2436)- 0.3848(0X0036) 

DXER24=-290. 9098+ 0.9820(0X2435)+ 1 . 84 2 2 ( D Y 24 35 ) 
+ 7 . 4308 ( XXLAT ) + 1 . 3908 ( DX1 224 ) 

- 20 . 7300 (VX0024)+ 1 2 . 9281 ( V Y1 224 ) 

- 14. 9940( VY0036) 

DYER24= 331.9491+ 5 . 4 1 0 3 ( XXL AT ) - 0 . 5 1 06 ( D Y 24 36 ) 

- 0 .6443( DX 1 2 24 ) - 4 . 7 0 7 8 ( V Y 00 1 2 ) 

- 1 . 1 65 8 ( XXLON ) 

DXER36=-431 . 3265+ 2 . 2 2 94 ( D Y 24 36 ) + 2 . 2 99 9 ( DX 1 2 24 ) 
+ 11 . 7478(XXLAT) + 2 1 . 66 34 ( '/ Y 1 2 24 ) 

+ 2.4042(DX2436)- 1.3375(0X0036) 

- 19.5337(VY0Q36) 

DY ER 36= 445.8439- 0.5419 ( 0X0024)- 2 . 76 2 1 ( DY 24 36 ) 
+ 5 . 6 707 ( XXLAT ) + 1 9 . 5 1 6 3 ( V Y 1 2 36 ) 

- 1 . 42 33 ( XXLON ) 



40 . 0% 

46.6% 

39.3% 

52.0% 

c ( 

to 

VARIANCE 

EXPLAINED 

50 . 7% 

51 .9% 

44 .2% 

31.4% 

3 7.8% 

42 .9% 
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B. FORWARD/BACKWARD INTEGRATION (1977-78 DATA) 



The regression equations which follow were produced by 



SPSS using predictors derived from forward and backward in- 
tegration of the TCM. The data sample consisted of 33 cases, 
or runs of theTCM in 1977 and 1978. Predictors in the equa- 
tions are listed in order of decreasing significance. 

1. TCM Integrated to 72 HR, and -36 HR % 

VARIANCE 

REGRESSION EQUATIONS " EXPLAINED 

D X ER 1 2 = 133.7858- 1 . 1 4 34 ( XXLAT ) +4 9 . 0 1 6 2 ( V X 1 2 36 ) 35.1% 

- 0.8544( JULDAY)- 1 . 3384 ( DX0036 ) 

+ 0.8434(DYM200)+ 0 . 2349 ( DXM6M2 ) 

+ 8 . 1 9 74 ( V YM5M4 ) - 1 . 9 5 3 2 ( V X6 0 7 2 ) 

+ 0.2739(DY5072) 

DYER12 = -323.6042- 16 . 76 34 ( V Y00 1 2 ) +1 1 . 781 2 (XXLAT) 86.5% 

+ 38.0056(VY0060)- 2 . 3 1 1 8 ( V Y 00 7 2 ) 

- 0.7205(DY2448)+ 2 1 . 82 8 2 ( V X 24 36 ) 

- 34 . 3595 ( VX0036 ) + 0 . 6 2 2 3 ( J UL DA Y ) 

- 0 . 3 5 7 9 ( DXM5M4 ) + 0 . 24 9 9 ( D XM4M 2 ) 

DXER24=1499 . 882+ 4 . 0720 ( VYM200 ) - 1 . 01 76 ( DYM6M4 ) 88.9% 

- 2 . 7329 ( JULDAY ) + 8 . 980 8 ( V X 1 2 3 5 ) 

+ O . 0 30 2 ( DXM6M2 ) - 4 . 7254 (XXLON ) 

+ 1 . 3 5 1 3 ( DY2436 )+ 1 9 . 7 88 1 ( V XM40 0 ) 

D YER24= 1 32. 901 5 + 1 3 . 1693 (XXLAT) - 0 . 46 89 ( DYOO 2 4 ) 7 7 . 3% 

- 12.0194(VX6072)- 1 . 2499 ( XXLON ) 

+ 10 . 8362 ( VX2436 ) - 0 . 3 1 6 1 ( DXM6M4 ) 

+ 12.5733(VY1224)- 0 . 5091 ( DY 3648 ) 

+ 0.5794(0X4372)- 1 2 . 1 87 9 ( V X00 36 ) 

DXER36 = 231 3. 2390 + 14. 1 0 2 5 ( V YM200 ) - 1 .4677 ( DYM5M4 ) 89.6% 

- 4.9703(JULDAY)+ 1 3 . 830 1 ( V X 1 2 36 ) 

+ 18.3607 (VXT14M2) - 4 . 884 1 ( XX LON ) 

- 12 . 2852 ( XXLAT) 

D Y ER 36= 208. 2644 + 13. 0 7 1 5 ( XXLAT ) - 2 . 1 86 3 ( X X LO N ) 74.9% 

- 2. 3350 ( VX4872) + 0.2046(VY2436) 

- 0. 7327 ( DXM5M4) + 12.3656 ( VYM6M4) 

+ 0.9038( JULDAY) 
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DXER 4 8=3766. 9010+ 3 . 96 98 ( V YM2Q0 ) - 7 . 5 0 5 3 ( J UL DA Y ) 
+ 2. 9278(VX4872)- 2 7 . 26 2 3 ( X XLA T ) 

- 10. 1 24 8 ( VY0048)- 18. 3859 ( VYM6M4 ) 

+ 24. 7539 ( VXM4M2 ) - 7. 241 7( XXLON) 

+ 27. 7 76 5(7X12 36 )- 1.3987( 0X3648) 

DY ER43= 587.8164- 1 . 2470 ( DXM6M2 ) - 1 3 . 4294 ( VY2436 ) 

- 5 . 9702 ( VX4872 ) + 91 . 4204 ( VYM424 ) 

- 4 . 34 1 8 ( XXLON ) + 0.2002(0X3648) 

+ 8 . 8435 ( XXLAT) + 1 . 8 326 ( J ULDA Y ) 

+ 2 . 02 70 ( DXM4M2 ) - 1 . 424 5 ( DYM2 1 2 ) 

+ 20 . 34 20 ( V YM6M4 ) - 1 . 3766 ( DYM200 ) 

DXER60=4803.057 +19.0639(VYM200)-10. 5436 (J ULDA Y) 
+ 8. 5565( VX4372)- 4 2 . 8 3 14 ( X X LA T ) 

- 46 . 1492(VY0024)- 6 . 9290 ( XXLON ) 

+ 2 . 3348 ( DXM4M2 ) - 1 . 66 72 ( DYM6M4 ) 

D Y ER60= 364.9735+ 0 . 6 7 1 1 ( D XM6 0 0 ) + 4 . 1 587 ( J UL DA Y ) 
+ 21 . 06 52 ( XXLAT ) + 7 5 . 96 30 ( V Y 00 2 4 ) 

+ 2.0039(DYM6M2)- 7 . 1 881 ( XXLON ) 

- 2 . 7226 ( DXM6M4 ) - 2 . 060 2 ( DYM2 1 2 ) 

DXER 72=3715.035 +87 . 3945 ( VYM200 ) - 6 . 6852 ( DYM6M4 ) 

- 1 2 . 76 0 3 ( JULDAY)- 115 .0140( VYC048) 

+ 42 . 7842 ( VXM6M2 ) - 39 . 16 36 ( XXLAT) 

- 0.8422(0X6072) 

DYER 72 = 1214. 7320-26. 6951 ( VYM200)- 7 . 6 1 89 ( XXLON ) 

- 2 . 792 1 ( DXM6M4 ) + 6 2 . 8756 ( VYM6M4 ) 

+ 38.9121 (VY00 24J+ 3 . 1 0 78 ( J U LDA Y ) 

+ 1 . 89-6 1 ( DY 3648 ) + 2 1 . 5 76 3 ( VX 3648 ) 

- 23.6336(VX0012) 



2 . TCM Integrated to 60 HR, and -36 HR 



REGRESSION EQUATIONS 

DXER 1 2= 16 7.5636- 1 . 1 6 79 ( XX LAT ) +4 2 . 0 954 ( V X 1 2 36 ) 

- 0 . 95 5 6 ( JULDAY)- 1 . 1 599 ( DX0036 ) 

+ 0 . 8736 ( DYM200 )+ 0 . 2185 ( DXM6M2 ) 

- 5 . 6 37 7 ( VYM6M4) 

DYER1 2 = - 344. 2635- 14. 9407 (VYOO 12 ) + 12. 0804 ( XXLAT) 

+ 26 . 5455 ( VY0060 ) + 1.6 737( 0X2436 ) 

- 26 . 5420(VX0036 )- 0 . 660 1 ( DY 3648 ) 

+ 0 . 6 4 7 3 ( JULDAY )- 0 . 3001 ( DXM6M4 ) 

+ 4 . 9 1 1 6 ( VXM 200 ) 



9 2.4% 



81.0% 



33.6% 



30.9% 



95 . 7% 



93.2% 



VARIANCE 

EXPLAINED 

32 . 2 % 



86 . 4 % 
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DXER24=1499.8820+ 4 . 0720 ( VYM200 ) - 1 . 01 76 ( 0YM6M4) 88.9% 

- 2. 7 32 9 ( JULDAY ) + 8.9808(7X1236 ) 

+ 0 . 0302 ( 0XM6M2 ) - 4 . 7254 ( XXLON ) 

+ 1 . 351 3 ( DY2436 ) + 19 . 7881 ( VXM400 ) 

DYER24 = - 56 .4323+12 . 9000( XXLAT) - 0 . 3 90 5 ( D Y 00 24 ) 72.5 % 

+ 9.0748(VY4860)- 0 . 3743 ( DXM6M4 ) 

+ 0. 5430( JULDAY)- 0 . 96 9 5 ( X X LON ) 

DXER36=2312. 2390+14. 1 02 5 ( V YM200 ) - 1 . 46 77 ( DYM6M4 ) 89.6% 

- 4.9703(JULDAY)+ 1 3 . 8 30 1 ( V X 1 2 36 ) 

+ 1 8 . 86 07 ( VXM4M2 ) - 4 . 884 1 ( XXLON ) 

- 1 2 . 285 2 ( X X LAT ) 

0 Y ER 36 = 518.1446 + 14. 5509 (XXLAT) - 3'. 906 7 ( XXLON ) 73.4% 

- 4 . 4077 ( VX4860 ) + 1 . 0088( JULDAY ) 

+ 1 . 2 5 97 ( DYM6M2 ) - 1 . 6 44 5 ( DXM6M4 ) 

+ 1 9 . 1 6 52 ( V XM6 00 ) + 9 . 26 26 ( V Y 1 224 ) 

- 1 3 . 881 3 ( VYM400) 

DXER48=3901 . 1060+ 0 . 1 4 32 ( V YM200 ) - 7 . 1 76 5 ( JULDAY ) 92.2% 

- 26 .0746 (XXLAT) + 33 . 2 885 ( V X 1 2 36 ) 

- 8 . 5896 ( XXLON ) + 0 . 489 3 ( D Y24 36 ) 

+ 2 . 0444 ( DXM4M2 ) - 1.3676(0X3648) 

- 1 . 3292( DYM6M4) 

DYER48=1228.8470- 3. 0296 ( DXM6M2) - 9 . 980 1 ( VY24 36 ) 81.4% 

+ 56.0325(VYM6M2)+ 1 .4210 ( JULDAY ) 

+ 50.1396 ( VY0024) - 4 7 . 0 803 ( V YM2 1 2 ) 

- 5 . 83 7 5 ( XXLON) + 2 . 5 5 1 8 ( 0 XM4M 2 ) 

+ 43.1560(VXM600)- 7 . 6 143 ( VX2436 ) 

DXER60 = 5564.0920- 7 . 204 1 ( VYM200 ) - 10 . 5 756 ( JULDAY ) 90.5% 

- 55.6857(XXLAT) - 1 247 . 4050 ( VY0024 ) 

+ 21 . 4406 ( VX4860 ) - 9 . 8929 ( XXLON ) 

+ 1833.842 ( V Y 00 36 ) + 2 . 3 39 2 ( DXM4M 2 ) 

—620 ,6515(VY2436 )+ 26.8308(VY2443) 

DYER50= 364.9736+ 0 . 6 7 1 1 ( DXM6 00 ) + 4 . 1 587 ( J U L DA Y ) 30.9% 

+ 21 . 06 52 ( XXLAT ) + 7 5 . 9 6 30 ( V Y 00 2 4 ) 

+ 2 . 0389 ( DYM6M2 ) - 7 . 1881 ( XXLON ) 

- 2 . 7226 ( DXM6M4 ) - 2 . 06 0 2 ( D YM2 1 2 ) 
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3 . TCM Integrated to 48 HR, and -36 HR 



<y 



VARIANCE 

REGRESSION EQUATIONS EXPLAINED 

D X ER 1 2 = 167.5636- 1. 16 78 ( XXLAT ) +4 2 . 0 95 4 ( V X 1 2 36 ) 32.2% 

- 0.9556( JULDAY)- 1 . 1 5 99 ( D XOO 36 ) 

+ 0.8736(DYM200)+ 0 . 2 1 84 ( DXM6M2 ) 

- 5 . 6 37 7 ( V YM6M4 ) 



DYER12 = -283 . 4232- 8 . 4050 ( V Y 001 2 ) + 9 . 2910 ( XXLAT ) 85.5% 

+ 12.9974(VY 1236 )+ 1 . 2 96 2 ( D X 2 4 36 ) 

- 1 7 . 6 26 7 ( V X00 36 ) + 0 . 6 3 5 5 ( J ULDA Y ) 

+ 5 . 6809 ( VYM6M4 ) - 0 . 2878 ( DXM6M4 ) 

DXER24=1499.8820+ 4.0720(VYM200)- 1 . 01 76 ( DYM6M4 ) 38.9% 

- 2 . 7 32 9 ( JULDAY)+ 8 . 9308 ( V X 1 2 36 ) 

+ 0 . 0302 ( DXM6M2 ) - 4 . 7254 ( XXLON ) 

+ 1 . 3 5 1 3 ( DY2436) + 19. 7881 (VXM400) 

DYER24= 129.5675+1 1 . 7946 (XXLAT) - 0 . 4 3 7 5 ( D Y 00 24 ) 75.9% 

+ 10.8043(VY1236)- 0 . 5 34 5 ( DXM6M4 ) 

+ 7 . 94 36 ( V YM6M4 ) - 1 . 4 7 70 ( XXLON ) 

+ 0.5032( JULDAY) 

DXER 36 = 2 31 3 . 2390 + 14. 1 0 2 5 ( V YM 200 ) - 1 . 46 77 ( D YM6M4 ) 89.6% 

- 4.9703(JULDAY ) + 1 3 . 83 0 1 ( V X 1 2 36 ) 

+ 18.8607(VXM4M2)- 4 . 834 1 ( XXLON ) 

- 1 2 . 285 2 ( XXLAT) 

D YER36= 130 .4497 + 1 3. 4500( XXLAT) - 1 . 93 1 9 ( X X LO N ) 74.5% 

- 0.8659(DXM6M4)+13.7236(VYM6M4) 

+ 0 . 91 30 ( JULDAY) 

DXER48=3901 . 1060+ 0 . 1 4 32 ( V YM 2 00 ) - 7 . 1 76 5 ( JU L D A Y ) 92.2% 

- 26 .0746 ( XXLAT ) +33.2885 (7X1236 ) 

- 8 . 5896 (XXLON ) + 0 . 489 3 ( D Y 24 36 ) 

+ 2 .0444( DXM4M2 )- 1.36 76 ( 0X3648) 

- 1 . 32 92 ( DYM6M4) 

DYE R4 8=1228. 847 - 3 . 0296 ( DXM6M 2 ) - 9 . 980 1 ( V Y 24 36 ) 81.4% 

+ 5 6 . 0 3 2 5 ( V Y M 5 M 2 ) + 1 . 42 10 ( JULDAY ) 

+ 50. 1396(VY0024)-47.0803(VYM212) 

- 5.8375(XXLQN) + 2 . 5 5 1 8 ( DXM4M2 ) 

+ 43. 1 5 6 0 ( VXM600)- 7.6143( 7X2436 ) 
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4 . TCM Integrated to 36 HR, and -36 HR 
REGRESSION EQUATIONS 

D X E R 1 2 = 167.5686- 1 . 16 79 ( XXLAT) +4 2 . 0 9 5 4 ( V X 1 2 36 ) 

- 0.9555 (JULDAY)- 1 . 1 5 99 ( DXOO 36 ) 

+ 0.8786(DYM200)+ 0 . 2185 ( DXM6M2) 

- 5 . 6 3 7 7 ( V YM5M4 ) 

DYER12 = -283 .4232- 8 . 40 4 9 ( V Y 00 1 2 ) + 9 . 29 1 0 ( XXLAT ) 

+ 12.99 74 (VY1236 )+ 1 . 2952 ( DX2436 ) 

- 17 .6267 (VX0036 )+ 0 . 6 355 ( JULDAY ) 

+ 5 . 6309 ( V YM6M4 ) - 0 . 2 8 78 ( DXM6M4 ) 

DXER24=1499.8820+ 4 . 0 720 ( V YM200 ) - 1 . 0176 ( DYM6M4) 

- 2 . 7329 ( JULDAY ) + 8 . 9 808 ( V X 1 2 36 ) 

+ 0 . 0 302 ( DXM6M2 ) - 4 . 7254 ( XXLON ) 

+ 1.3513(DY2436)+19.7881(VXM400) 

DY ER24= 129. 5675 + 1 1. 7946 (XXLAT) - 0 . 4375 ( DY0024 ) 
+ 10 .8043(VY1236 ) - 0 . 5 34 5 ( BXM6M4 ) 

+ 7.9436(VYM6M4)- 1 . 5 770 ( XXLON ) 

+ 0.5032( JULDAY) 

DXER 36 =2 3 13. 2 3 90+1 4. 1025 ( VYM200)- 1 . 46 7 7 ( DYM6M4 ) 

- 4.9703(JULDAY)+13.8301(VX1236) 

+ 18.8607 ( VXM4M2) - 4 . 834 1 ( X XLON ) 

- 12 . 2852 ( XXLAT ) 

DYER36= 130. 4497 + 13. 4500 ( XXLAT) - 1 . 9 3 1 9 ( XX LO N ) 

- 0 . 86 5 9 ( DXM6M4 ) + 1 3 . 7 2 36 ( VYM6M4 ) 

+ 0 . 9 1 30 ( JULDAY) 



% 

VARIANCE 

EXPLAINED 

82 . 2% 



35.5% 



88 . 9 % 



75.9% 



89.6% 



74 .5% 
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